#A long-term restoration experiment at Rancho Marino Reserve.
To explore existing evidence and data empirically and assess future research activities including sampling for efficicacy of intervention on recruitment of natives into coastal grassland with extensive exotic plant representation.
The control macroplots were not used.
From data sheet columns: Plot 12x12m(12 setup & originally sampled, but as you noted control sites weren’t uses) Subplot 3x3m (4 per plot)
Quadrat 1x1m (5 per Subplot)
Don Canestro In summary, 4 years of sampling 1 before and three after planning. I also have data on plug grow out and which spp. seeds had best germination rates and plot images from 2006 and 2009.
Between 2002 before data and 2006 data collection was kill time. I was told by restoration experts 3 years of death on a plot was good to remove seed bank and reduce competition
Below are descriptions of the attached files. “harding grassexpt3_05.doc” describes the experiment “HARDINGGRASSPLOTS.xls” illustrates the plot layout. “Percent cover data2002.xls” is before treatment data Other files self evident.
library(tidyverse)
library(plotly)
#2009 only
data <- read_csv("data/2009.data.csv")
data
## # A tibble: 160 × 73
## data plot plot.treatment subplot subplot.treatment.number
## <int> <int> <chr> <chr> <int>
## 1 2009 1 plastic D 1
## 2 2009 1 plastic D 1
## 3 2009 1 plastic D 1
## 4 2009 1 plastic D 1
## 5 2009 1 plastic D 1
## 6 2009 2 till.spray A 1
## 7 2009 2 till.spray A 1
## 8 2009 2 till.spray A 1
## 9 2009 2 till.spray A 1
## 10 2009 2 till.spray A 1
## # ... with 150 more rows, and 68 more variables: subplot.treatment <chr>,
## # quadrat <int>, max.height <int>, tallest.species <chr>, `small mammal
## # holes` <int>, bare <int>, thatch <int>, Aica <int>, Anar <int>,
## # Avena <int>, Baca <int>, Bado <int>, Bapi <int>, Brca <int>,
## # Brdi <int>, Brho <int>, Brmi <int>, Care <int>, Carp <int>,
## # Cheno <int>, Civu <int>, Coca <int>, Con <int>, Cyper <int>,
## # Cyper0sp2 <int>, Daca <int>, Elgl <int>, Eltr <int>, flat0sedge <int>,
## # Gedi <int>, Ger0sp <int>, Gna <int>, Hepu <int>, Hobr <int>,
## # Hypo <int>, Isce <int>, Jubu <int>, Jupa <int>, Junc.sp <int>,
## # Junc.sp.2 <int>, Juoc <int>, Lomu <int>, Lope <int>, Lyhy <int>,
## # Malva <int>, Mepo <int>, Myopo <int>, Napu <int>, Phaq <int>,
## # Poly <int>, Rasa <int>, Ruac <int>, Rucr <int>, Sil <int>, Sima <int>,
## # Sonch <int>, Spma <int>, Trifol0sp <int>, forb <int>, grass <int>,
## # Vela <int>, Ver <int>, Vumy <int>, total.cover <int>,
## # density.Daca <int>, density.Elgl <int>, density.Napu <int>,
## # total.native.density <int>
#str(data)
#totals
p <- ggplot (data, aes(plot.treatment, total.cover)) + geom_boxplot() + facet_wrap(~subplot.treatment)
ggplotly(p)
p <- ggplot (data, aes(plot.treatment, total.native.density)) + geom_boxplot() + facet_wrap(~subplot.treatment)
ggplotly(p)
p <- ggplot (data, aes(plot.treatment, total.native.density, color = subplot.treatment)) + geom_boxplot() + geom_point()
ggplotly(p)
#for each native
p <- ggplot (data, aes(plot.treatment, density.Daca)) + geom_boxplot() + facet_wrap(~subplot.treatment)
ggplotly(p)
p <- ggplot (data, aes(plot.treatment, density.Elgl)) + geom_boxplot() + facet_wrap(~subplot.treatment)
ggplotly(p)
p <- ggplot (data, aes(plot.treatment, density.Napu)) + geom_boxplot() + facet_wrap(~subplot.treatment)
ggplotly(p)